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State variables, health parameters, and actuators 
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FIG. 2 


Sens ed outputs and standard deviation as percent of operating point trim va lues 
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Estimated auxiliary parameters 
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FIG. 4 


Auxiliary parameter squared estimation errors 
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Health parameter % squared estimation errors (reduced noise) 
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OPTIMIZED TUNER SELECTION FOR 
ENGINE PERFORMANCE ESTIMATION 

CROSS-REFERENCE TO RELATED 

APPLICATION 5 

This application claims benefit from U.S. Provisional 
Patent Application No. 61 / 247 , 325 , entitled “OPTIMAL 
TUNER SELECTION,” filed on Sep. 30 , 2009 , which is 
hereby incorporated by reference in its entirety. 

BACKGROUND 

An emerging approach in the field of aircraft engine con- 
trols and health management is the inclusion of real-time 
on-board models for the in-flight estimation of engine perfor- 
mance variations. This technology, typically based on Kal- 
man filter concepts, enables the estimation of unmeasured 
engine performance parameters that can be directly utilized 2 o 
by controls, prognostics and health management applica- 
tions. A challenge which complicates this practice is the fact 
that an aircraft engine’ s performance is affected by its level of 
degradation, generally described in terms of unmeasurable 
health parameters such as efficiencies and flow capacities 25 
related to each major engine module. Through Kalman filter- 
based estimation techniques, the level of engine performance 
degradation can be estimated, given that there are at least as 
many sensors as parameters to be estimated. However, in an 
aircraft engine the number of sensors available is typically 30 
less than the number of health parameters presenting an 
under-determined estimation problem. A common approach 
to address this shortcoming is to estimate a sub-set of the 
health parameters, referred to as model tuning parameters. 
While this approach enables on-line Kalman filter-based esti- 
mation, it can result in “smearing” the effects of unestimated 
health parameters onto those which are estimated, and in turn 
introduce error in the accuracy of overall model-based per- 
formance estimation applications. 40 

Recently, a new method has been presented based on sin- 
gular value decomposition that selects a model tuning param- 
eter vector of low-enough dimension to be estimated by a 
Kalman filter. The model tuning parameter vector, defined as 
q, was constructed as a linear combination of all health 45 
parameters, h, given by the equation 

q=V*h, (1) 

where the transformation matrix, V*, is selected applying 
singular value decomposition to capture the overall effect of 50 
the larger set of health parameters on the engine variables as 
closely as possible in the least squares sense. 

SUMMARY 

55 

A new linear point design technique which applies a sys- 
tematic approach to optimal tuning parameter selection is 
presented. This technique defines a transformation matrix, 
V*, used to construct a tuning parameter vector which is a 
linear combination of all health parameters, and of low 60 
enough dimension to enable Kalman filter estimation. The 
new approach optimally selects the transformation matrix, 
V*, to minimize the theoretical steady-state estimation error 
in the engine performance parameters of interest. There is no 
known closed form solution for optimally selecting V* to 65 
satisfy this objective. Therefore, a multivariable iterative 
search routine is applied to perform this function. 


2 

BRIEF DESCRIPTION OF THE SEVERAL 
VIEWS OF THE DRAWINGS 

Objects and advantages together with the operation of the 
invention may be better understood by reference to the 
detailed description taken in connection with the following 
illustrations, wherein: 

FIG. 1 illustrates a flow chart for performing an iterative 
optimal search; 

FIG. 2 illustrates a table of state variables, health param- 
eters, and actuators; 

FIG. 3 illustrates a table of sensed outputs and standard 
deviation as percent of operating point trim values; 

FIG. 4 illustrates a table of estimated auxiliary parameters; 

FIG. 5 illustrates a table of auxiliary parameter squared 
estimation errors; 

FIG. 6 illustrates a first graph of tuner comparisons; 

FIG. 7 illustrates a second graph of tuner comparisons; 

FIG. 8 illustrates a third graph of tuner comparisons; 

FIG. 9 illustrates a fourth graph of tuner comparisons; 

FIG. 10 illustrates a table of health parameter squared 
estimation errors at nominal noise levels; 

FIG. 11 illustrates a table of health parameter squared 
estimation errors at reduced noise levels; and 

FIG. 12 illustrates a graph of tuner impact on estimator 
response. 

DETAILED DESCRIPTION 

Reference will now be made in detail to exemplary 
embodiments of the present invention, examples of which are 
illustrated in the accompanying drawings. It is to be under- 
stood that other embodiments may be utilized and changes 
may be made without departing from the respective scope of 
the present invention. 

The following includes definitions of selected terms 
employed herein. The definitions include various examples 
and/or forms of components that fall within the scope of a 
term and that may be used for implementation. The examples 
are not intended to be limiting. Both singular and plural forms 
of terms may be within the definitions. 


V A xh ,A xq 

B > B x/»> B xg 

c,c x „c^ 
D,F,F*,F^ 
G, L, M, N 

System Matrices 

C-MAPSS 

Commercial Modular Aero-Propulsion 
System Simulation 

Fn 

Net Thrust 

G x, G a , G, 

Estimation bias matrices 

H 

Matrix which relates health parameter 
effects to steady-state engine outputs 

HPC 

High pressure compressor 

HPT 

High pressure turbine 

I 

Identity matrix 

Koo 

Kalman filter gain 

LPC 

Low pressure compressor 

LPT 

Low pressure turbine 

MAP 

Maximum a posteriori 

Nf 

Fan speed 

Nc 

Core speed 


Health & auxiliary parameter covariance matrices 

P24 

HPC inlet total pressure 

Ps30 

HPC exit static pressure 

P P P 

xh,k 5 ’ 

Covariance matrices of estimated parameters 

Poo 

Kalman filter state estimation covariance matrix 

Q> Q xhi Q xq 

Process noise covariance matrices 

R 

Measurement noise covariance matrix 
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4 

-continued 


SmLPC 

T24 

T30 

T40 

T48 

T50 

V* 

VSV 

VBV 

Wf 

w z 

h 

q* 

% 

v* 

W k> W h,k> W xh,k 
x k 

X xh,k 

X xq,k 

y k 

Zk 

F'xqjc 

Subscripts 


LPC Stall margin 

HPC inlet total temperature 

HPC exit total temperature 

Combustor exit temperature 5 

Exhaust gas temperature 
LPT exit temperature 
Transformation matrix relating h*. to q* 

Variable stator vane 
Variable bleed valve 

Fuel flow 10 

Auxiliary parameter weighting matrix 

Health parameter vector 

Kalman filter tuning parameter vector 

Actuator command vector 

Measurement noise vector 

Process noise vectors 1 ^ 

State vector 

Augmented state vector (xk and hk) 

Reduced order state vector (xk and qk) 

Vector of measured outputs 

Vector of unmeasured (auxiliary) outputs 

residual vector (estimate minus its expected value) 


k 

xh 

xq 

ss 

Superscripts 


Discrete time step index 
Augmented state vector (x and h) 
Reduced order state vector (x and q) 
Steady-state value 


25 


f Pseudo- inverse 

Estimated value 

~ Error value 

Mean value 

T Transpose 30 

Operators 


EH 
tr {'} 

SSEE(-) 

WSSEE(-) 

IiIf 


Expected value of argument 
Trace of matrix 

Sum of squared estimation errors 
Weighted sum of squared estimation errors 
Matrix Frobenius norm 


The discrete linear time-invariant engine state space equa- 
tions about a linear design point are given as: 


40 


=Ax k +Bu k +Lh k +w k ; 


y k = Cx k +Du k +Mh k +v k 


z k =Fx k +Fu k +Nh k (2) 

where k is the time index, x is the vector of state variables, u 
is the vector of control inputs, y is the vector of measured 
outputs, and z is the vector of auxiliary (unmeasured) model 
outputs. The vector h represents the engine health parameters, 
which induce shifts in other variables as the health parameters 
deviate from their nominal values. The vectors w and v are 
uncorrelated zero-mean white noise input sequences. Q will 
be used to denote the covariance of w, and R to denote the 
covariance of v. The matrices A, B, C, D, F, G, L, M, andN are 
of appropriate dimension. The health parameters, represented 
by the vector h, are unknown inputs to the system. They may 
be treated as a set of biases, and are thus modeled without 
dynamics. With this interpretation Eq. (2) can be written as: 


r x k 

yt = [ c M 1 , +Du k +V k 

TAJ 

x xh,k 

= C x h x xh, k + Du k + v k 

U = [F ivTM + Gn* 

~zr [ -M 

x xh,k 

= F xh x xh , k + Gu k 


The vector w x/2 is zero-mean white noise associated with 
the augmented state vector, [x 7 h 7 | r , with a covariance of Q xh . 
W xh consists of the original state process noise, w, concat- 
enated with the process noise associated with the health 
parameter vector, w h . 


w xhjc = 



(4) 


The eigenvalues of A xh consist of the original eigenvalues 
of A plus an additional dim(h) eigenvalues located at 1 .0 on 
the unit circle due to the augmentation. Thus, the new aug- 
mented system given in Eq. (3) has at least as many eigenval- 
ues located on the unit circle as there are elements of h. Once 
the h vector is appended to the state vector, it may be directly 
estimated, provided that the realization in Eq. (3) is observ- 
able. Using this formulation, the number of health parameters 
that can be estimated is limited to the number of sensors, the 
dimension of y. Since in an aircraft gas turbine engine there 
are usually fewer sensors than health parameters, the problem 
becomes one of choosing the best set of tuners for the appli- 
cation. A methodology is presented for the optimal selection 
of a model tuning parameter vector, q, of low-enough dimen- 
sion to be estimated by a Kalman filter, while minimizing the 
estimation error in the model variables of interest. The steps 
in this process include construction of the reduced-order state 
space model, formulation of the Kalman filter estimator, cal- 
culation of the mean sum of squared estimation errors, and 
optimal selection of the transformation matrix to minimize 
the estimation error. 

The first step is to construct a reduced-order state space 
model. The model tuning parameter vector, q, is constructed 
as a linear combination of all health parameters, h, given by 

q=V*h (5) 

where qEP m , hEP^, m<p, and V* is an mxp transformation 
matrix of rank m, applied to construct the tuning parameter 
vector. An approximation of the health parameter vector, h, 
can be obtained as 


h=V*\ (6) 

where V*^ is the pseudo -inverse of V*. Substituting Eq. (6) 
into Eq. (3) yields the following reduced order state space 
equations which will be used to formulate the Kalman filter: 


'*a+i 
. h+i 


A L 
0 1 


A xh x xh,k 



w xh,k 


(3) 


**+l 

A 

urn 

r x k r b i 

Wk 

. #t+l J 

.0 

/ 

+ n u k + 

L qk L o J 



A xq x xq,k B xq w xqk 


— A X qX X q k 4 - B X qll k + W X q, k 


( 7 ) 


= A x hX x h,k + B xh u k + w xfhk 
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-continued 

y k = [C_MV^]^‘ ] + Du k + v k 

C xq 

x xq,k 

— Cxq x xq,k + Du k + V k 


Z k = [F NV*]\ Xt ] + Gu k 

~— F iiii 

t xq 

x xq,k 10 

= F X q x xq,k + Gu k 


6 

ance. The estimation error bias vectors are equivalent to the 
mean estimation error vectors defined as: 

Xxhjc = E[x xhJl ] = E\x xhtk - x xKk ] ( 1 4 ) 

Ik =E[zk] = E[zk - Zk ] 

where the operator E[*] represents the expected value of the 
argument. The variance of the estimates can be found by 
constructing their respective estimation covariance matrices: 


The state process noise, w xq , and its associated covariance, 
Q xr/ , for the reduced order system are calculated as: 


E — E[(x x h, k — E\x x h jk \){x x h jk — E\x xk k ]) ] 

P u = E[(Zk -E[zk]){z k -E[zic]f] 


(15) 


= 


Qxq = 


I 

0 ■ 

/ 

0 ■ 

■ Wk 

0 

y* W xh,k = 

.0 

v* _ 

. w h,k . 


I 0 

0 v* 


Qxh 


/ 0 
0 V* 


(8) 


Next, the Kalman filter estimator is formulated. Here, 
steady-state Kalman filtering may be applied. This means that 
while the Kalman filter is a dynamic system, the state estima- 
tion error covariance matrix and the Kalman gain matrix are 
invariant — instead of updating these matrices each time step 
they are held constant. Given the reduced order linear state 
space equations shown in Eq. (7), the state estimation error 
covariance matrix, P rjo , is calculated by solving the following 
Ricatti equation: 

p^=a^a x /-a x ^c x /{c x ^c x / + ry 

i C xq P«A xq T +Q xq (9) 

The steady- state Kalman filter gain, K TO , can then be calcu- 
lated as follows: 

=P O0 C xq T (C xq P a:> C xq T +R )- 1 (10) 


Diagonal elements of the covariance matrices will reflect the 
20 variance in individual parameter estimates, while off-diago- 
nal elements reflect the covariance between parameter esti- 
mates. The overall sum of squared estimation errors (SSEE) 
can be obtained by combining the estimation error bias and 
estimation variance information as 

25 S S EE{£ xh j)=\ h /X xKk +tr{Pz kk } 

SSEE(f k )=2 k %+tr{P z%k } (16) 

where tr{*} represents the trace (sum of the diagonal ele- 
ments) of the matrix. As described herein, theoretical values 
30 for each error component will be derived assuming steady- 
state, open-loop (u=0) operating conditions. First, the estima- 
tion error bias is derived, followed by a derivation of the 
estimation variance. 

The estimation error biases, x xhk and can be analytically 

35 derived for an arbitrary health parameter vector, h, at steady - 
state operating conditions. This is done taking advantage of 
the following expected value properties at steady-state open- 
loop operating conditions 


and, assuming steady-state, open-loop operation (u=0), the 40 E[x k+ J -E[xjj-x ss 

Kalman filter estimator takes the following form 

E[hjj=h 


X xq ,k=^ x ^ xq ,k^i+Kjy k - C xq A x ^ xq ^_ x ) (1 1) 

The reduced order state vector estimate, x xq>k , produced by 
Eq. (11) can be used to produce an estimate of the augmented 45 
state vector, and the auxiliary parameter vector as follows: 


E[x xhJ J=x xh ^ 


E[y k J=y i 

E[Z]J=z s 


X%h,k 


/ 0 

o r f 


X xq.k 


Zk =[E Nr^]x xq , k 


( 12 ) 

50 


EM= 0 
Efw k ]=0 
E[v k ]=0 


The estimation errors in x xh ^ k and z k are defined as the 55 
difference between the estimated and actual values: 

X xh,k =X xh,k-X x hJc 

h=Zk-Zk (13) 60 

Due to the under-determined nature of the estimation prob- 
lem, it will be impossible for the Kalman filter estimator to 
completely restore all information when transforming q into 
h. As such, the Kalman filter will be a biased estimator, 
meaning the expected values of x hJc , and z k will be non-zero. 65 
The estimation errors can be considered to consist of two 
components: an estimation error bias, and an estimation vari- 


E[x xq ^=E[x xqk _J=I xq ^ 

E[x xKk j=\ Kss 

EffjJ=f ss (17) 

where the subscript “ss” denotes steady-state operation. By 
taking expected values of Eq. (2), x^, y^ and z ss can be 
written as functions of the health parameter vector h 

E[x k+ J =A E[x/J+B -EfiijJ+L -EfhjJ+Ef w k ] 
x ss =Ax ss +Lh 
x s =(I-A)~ l Lh 


( 18 ) 
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EfyJ= OEfxJ+D-EfuJ+M-EpiJ+EM 


y ss =Cx ss +Mh 


y s =(C(I-A)~ l L+M)h 


-continued 

[F NV^d-A^ + K^A^r 1 ...x 

t» = K ao [C(I-A)- i L+M]-- 

[F(I -A)~ l L + N] 


E[z k j=F-E[x k J+G-E[u k J+N-E[h k \ 


z s =(F(I-A)- l L+N)h (20) 

Next, by taking expected values of both sides of Eq. (1 1), the 
expected value of x xqk can be obtained as a function of y ss : 

E[x xq ^}=A xq E[x k _ i], . . 


The estimation error bias equations given in Eqs. (23) and 
(24) are functions of an arbitrary health parameter vector, h. 
As such they are representative of the parameter estimation 
error biases in a single engine, at a given point in its lifetime 
of use where its deterioration is represented by the health 
parameter vector h. The average sum of squared estimation 
error biases across a fleet of engines can be calculated as 


+K^Efy k j-C xq A xq -Efx xqJc _ l ]) 

X xg^s~^xq X xq,ss^'E <X) (y ss — C xq A xq X xqss ) 20 

£ xq ss=(I-A xq +K^ C^)- 1 K„y ss (21) 

Then, making the substitution y^=(C(I-A) _1 L+M)^ given in 
Eq. (19), the expected steady-state value of x xq ^ can be writ- 25 
ten as a function of h 

^HI-Axq+K- C xq A xq )~ l KUC(I-A)- 1 L+M)h (22) 

The steady-state augmented state estimation error bias can 
then be found, and partitioned into error bias information for 30 
the original state vector, x ss , and the health parameter vector, 

E^, by combining Eqs. (12), (14), (18) and (22) to yield 


X xh, fleet ~ — 

= E[tr{G xh hh T G T xh }} 
= trjG xh - E[hh T ] G T xh 


-- tr{G xh P h G T xh } 


l fleet = 


= E[tr{G z hh G z }] 
--tr^G z - E[hh T ] G T 7 ^ 


X xh,ss - E[x xh k - X xhk ] 


X xh,ss — X xh,ss X xh,s. 


(/ - A xq + K^C^A^y 1 


Koo [C(/ - A) -1 L+ M] ... - 


(I-A)- l L] 


(23) where the matrix P A , defined as E[hh r ], reflects a priori or 
35 historical knowledge of the covariance in the health param- 
eters across all engines. If available, it can be used to predict 
the sum of squared estimation errors biases as shown in (25) 
and (26). 

4Q Next, derivations are presented for the augmented state 
estimate and auxiliary parameter estimate covariance matri- 
ces, and P f?A: , respectively. These matrices will be cal- 
culated as a function of the reduced-order state vector esti- 
mation covariance matrix, P^^, which is defined as 


Exq, k ~ E\ \ x xq,k E\x xq , k \)[X X q t k E\x xq , k ]) 


where the vector e. xq , k is defined as the residual between x xq ^ 
at time k and its expected value. Since E[x„ 7>t | =x x? ^, e xqX 
can be obtained by subtracting Eq. (21) from Eq. (11) 


The steady-state auxiliary parameter estimation error bias can 
also be derived by combining Eqs. (12), (14), (20) and (22) to 
yield 


^xq,k — X xq,k I . X xqjk J 


: ^sq x xq , k - 1 3" E 00 (y k E X q^xq x xq,k - 1 ) ■ ■ ■ 


i = E[h ~ Zk ] 


{/^xq x xq,ss 3" ^oo C xq A xq X xq ss 


, = [F Nr f ]x xq , ss -Zss 


: (y^xq C xq A xq ) {x xq , k _ I X xq^s) E oa (y k y ss ) 


65 
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Making the substitutions and v*=y*- 

y„ yields 

^ q jc^ xq -K^ xq A- xq )^-i+K^v k (29) 

The estimation covariance matrix is then calculated as 

p m = E ^* e i,^ (30) 

— \Axq ~ KooC X qA X q]E[£ X q,i i -i£ X qj { _Y\ ... X 
[&xq - Koo C xq A xq ] T ... + 

[&xq ~ K 00 C xq A xq ]E[£ xqtk -. l vl]Kl 0 ... + 


Once Eqs. (25), (26), (33) and (35) are obtained, they may 
be used to analytically calculate the mean sum of squared 
estimation errors over all engines by combining the respec- 
tive estimation error bias and estimation variance information 
5 as previously shown in Eq. (16). The mean augmented state 
vector sum of squared estimation errors, SSEE(x xAi/fee ,), and 
the mean auxiliary parameter vector sum of squared estima- 
tion errors, SSEE(x xAj/?ee? ), become 

to 

SSEE(x* fl '«) = t h 'j leet +tr{P Ak } (36) 

= tr{G xh P h G T xh + P kkk } 


Koo E[v k £ xqji - 1 1 \-^xq C X qA X q\ ... 4 - 

K„E[v k vl)Kl 


SSEE& fleet) = ffleet + tr { P Xk) 

= i ’r{G z P h G[ + P u ) 


The substitutions 'E,[z xq<k _ l e xqJc _ t T ] =P^ > *_ j , and E[v t v/]=R 
can be made in the above equation. Since e ^ and v k T are 20 
uncorrelated, the substitution E[6 x? i _ 1 v i 7 ']=E[v*6^_ 1 7 ']=0 
can also be made, producing 

P xqJc^f-d-xq-Kcn C xq A xq JP %2J c -l[A xq -K <x> C xq A xq ] T + 

KJLKj (31) 

.25 

At steady-state operating conditions ^xqjc-G^xqjc Making 
this substitution in (31) produces the following Ricatti equa- 
tion which can be solved for P^^: 

KJiKj (32) JU 

It should be noted that V^ k obtained by solving (32) will be 
identical to P^produced via Eq. (9) if the system’ s actual state 
process noise covariance is identical to the assumed in the 

design of the Kalman filter. However, Q is often treated as a 35 
Kalman filter design parameter to provide acceptable 
dynamic response. For the purpose of this derivation, we have 
assumed a steady- state operating condition where the state 
variables and health parameters are invariant, and thus the 
actual system process noise is zero (i.e., w xA ^=0). In this case 40 
^xq,k w iH n °f equal P^. Once P^^ is obtained, it can be used 
to calculate P^, the covariance of x xh ^ k , which is defined as 


p - = 

/ 

0 

P m 

/ 

0 ■ 

xh,k 

0 

v* + 


0 



The augmented state vector estimation covariance given in 
Eq. (33) can be partitioned into covariance information for the 
original state vector, P^ (upper left comer of the P xh,k 
matrix), and the health parameter vector, Pfi ,k (lower right 
corner of the k matrix) 




55 


(34) 


The ¥ X q k matrix from Eq. (32) can also be used to calculate 60 
P f *, the covariance in the estimation of z k , which is equivalent 
toE[(VE[zJ)(VE[zJ) r ] 

PzjrtF NV*^] T (35) 

The variance in the estimates x xh>k and z k can be obtained from 65 
the diagonals of the covariance matrices produced by (33) and 
(35) respectively. 


If required, a weighted sum approach can be applied to 
normalize the contributions of individual auxiliary parameter 
estimation errors. This is often necessary as there may be 
several orders of magnitude difference between the auxiliary 
parameters of interest. A weighted sum approach prevents 
domination by individual parameters. In this study a diagonal 
auxiliary parameter weighting matrix, W z , is applied based on 
the inverse of auxiliary parameter variance (obtained from the 
main diagonal of the auxiliary parameter covariance matrix, 

PJ 


P z = [F(I -A)~ l L + N]P h [F(I - A)~ l L + N] T 


W z = 


Pz. 11 0 o r 1 

0 ■■. 0 

o o P^n . 


(37) 


W z is then applied to calculate of a “weighted” sum of auxil- 
iary parameter squared estimation errors given as 

WSSEE {f fleet )=tr{ WJG^GZ+P^] } (3 8) 

From Eqs. (23), (24), (33) and (35) it can be observed that 
both bias and variance are affected by the selection of the 
transformation matrix, V*. The sum of squared estimation 
error terms derived in this section give rise to an optimization 
problem: selecting V* to minimize the squared estimation 
error in the Kalman filter produced parameter estimates. This 
could include health parameter estimates, auxiliary param- 
eter estimates, or a combination of parameters. Although 
there is no known closed form solution for optimally selecting 
the V* matrix to satisfy the objective of minimizing estima- 
tion errors, a multi-parameter iterative search method has 
been developed to perform this task, and will be described in 
the next section. 

Prior to initiating the search for an optimal V*, specific 
system design information must be defined or obtained. This 
includes specifying the auxiliary parameters to be estimated, 
generating system state space equations at a fleet average 
(50% deteriorated) engine trim point, defining measurement 
noise covariance matrix, R, defining augmented state process 
noise covariance matrix, Q xh , and defining fleet average 
health parameter covariance, V h . 

Some additional clarification is provided regarding the 
selection of P A and Q xh as the distinction between these two 
ovariance matrices may not be immediately obvious. P A 
defines the expected health parameter covariance across all 
engines. It may be based on past knowledge gained from 
engine gas path analysis programs and/or historical studies of 
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engine module performance deterioration. Conversely, Q* A 
defines the expected process noise covariance in the state 
variables and health parameters of an individual engine, at a 
single discrete time step, k. The selection of Q Jch will directly 
impact the dynamic response and the variance of the esti- 
mates generated. 

After the necessary system information has been obtained, 
the search for an optimal transformation matrix to minimize 
the Kalman filter sum of squared estimation errors can com- 
mence. This is performed using the lsqnonlin function of the 
Matlab® Optimization Toolbox. This function applies an 
iterative search to find the least squares solution of a user- 
specified multivariable optimization problem. A flow chart 
depicting the steps in this optimal iterative search is shown in 
FIG. 1, and a further description of each step is given below. 

Upon startup, an initial random guess ofV* is generated. It 
is selected such that the matrix Frobenius norm ||V*|[p=l. This 
requirement is applied to help prevent convergence to a 
poorly scaled solution. The reduced order state-space model 
is then constructed, (Eq. (7)). The Kalman filter is then for- 
mulated by first calculating the estimation covariance matrix, 
(Eq. (9)), then calculating the Kalman gain matrix, K^, 
(Eq. (10)). Next, the sum of squared estimation errors (Eq. 
(36)), or weighted sum of squared estimation errors (Eq. (38)) 
is calculated. On each iteration the change in SSEE (or 
WSSEE) relative to the previous iteration is assessed to deter- 
mine if convergence within a user specified tolerance has 
been achieved. If convergence is not achieved then, V* is 
updated via the Matlab® lsqnonlin function, again requiring 
that ||V*|[p=l, and then returning to again construct the 
reduced order state- space equations. If convergence is 
achieved then he optimization routine returns the optimal 
value ofV*, and ends. 

The transformation matrix returned by the optimization 
routine may not be unique. Thus, different matrices can be 
found which produce a global minimum of the objective 
function. Experience has also shown that the optimization 
routine will usually return a V* matrix which satisfies, or 
nearly satisfies (i.e. within 5%), the global minimum of the 
objective function. However, in order to guard against poten- 
tial convergence to a local minimum, it is prudent to run the 
optimization routine multiple times, each time starting with a 
different initial guess for V*. This is only to assure the 
designer that the global minimum is achieved, not to produce 
a consistent V*. It should be emphasized that the optimal 
search for V* is only conducted off-line during the estimator 
design process . This calculation is not conducted as part of the 
on-line real-time Kalman filter implementation, and thus 
places no additional computational burden upon it. 

In an applied example of the methods described herein, a 
linearized cruise operating point extracted from the NASA 
Commercial Modular Aero -Propulsion System Simulation 
(C-MAPSS) high-bypass turbofan engine model is used to 
evaluate the new systematic tuner selection methodology. 
The linear model has two state variables, ten health param- 
eters, and three control inputs, all shown in FIG. 2. The 
model’s seven sensed outputs, and corresponding sensor 
noise standard deviation, are shown in FIG. 3. The auxiliary 
output parameters of interest to be estimated are shown in 
FIG. 4 . The linear model is used as the truth model for this 
application example. The model is run open-loop, so all con- 
trol inputs remain at 0, i.e., they do not deviate from the trim 
value for the linear model and no actuator bias is present. 
Deviations in all ten health parameters are assumed to be 
uncorrelated, and randomly shifted from their trim conditions 
with a standard deviation of ±0.02 (±2%). Since a parameter’ s 
variance is equal to its standard deviation squared, the health 
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parameter covariance matrix, P^, is defined as a diagonal 
matrix with all diagonal elements equal to 0.0004. 

Next, the estimation accuracy of the systematic approach 
for selecting Kalman filter tuning parameters will be com- 
5 pared to the conventional approach of selecting a sub- set of 
health parameters to serve as tuners (the seven health param- 
eters denoted with in FIG. 2), and the singular value 
decomposition approach to tuner selection. FIG. 5 shows a 
comparison of the theoretically predicted estimation errors 
10 (squared bias, variance, and total squared error) and experi- 
mentally obtained squared estimation errors for each of the 
three tuner selection approaches. T40 and T50 estimation 
errors are shown in squared degrees Rankine, and Fn and 
SmLPC estimation errors are shown in squared percent net 
15 thrust and squared percent stall margin respectively. The 
experimental results were obtained through a Monte Carlo 
simulation analysis where the health parameters varied over a 
random distribution in accordance with the covariance 
matrix, P A . The test cases were concatenated to produce a 
20 single time history input which was provided to the 
C-MAPSS linear discrete state space model given in Eq. (2), 
with an update rate of 15 ms. Each individual health param- 
eter test case lasted 30 s. 

At the completion of each 30 s test case, the health param- 
25 eter vector input instantaneously transitioned to the next test 
case. A total of 375 30 s test cases were evaluated, resulting in 
an 1 1,250 s input time hi story. Three separate Kalman filters 
were implemented using the three tuner selection approaches . 
The experimental estimation errors were determined by cal- 
30 culating the mean squared error between estimated and actual 
values during the last 10 s of each 30 s test case. The error 
calculation is based on only the last 1 0 s so that engine model 
outputs and Kalman estimator outputs have reached a quasi - 
steady-state operating condition prior to calculating the error. 
35 This ensures that the experimental results are consistent with 
the theoretically predicted estimation errors which were 
derived assuming steady-state operation. 

From FIG. 5 it can be seen that the theoretically predicted 
and the experimentally obtained squared estimation errors 
40 exhibit good agreement. If the number of random test cases 
were increased to a suitably large number, it is expected that 
the theoretical and experimental results would be identical. It 
can also be seen that all three estimators are able to produce 
unbiased estimates of the combustor exit temperature, T40; 
45 however, their estimates of LPT exit temperature, T50, net 
thrust, Fn, and LPC stall margin, SmLPC, are biased. The 
encouraging finding is that the new systematic approach to 
tuner selection significantly reduces the overall mean squared 
estimation error compared to the other two approaches. Rela- 
50 five to the conventional approach of tuner selection the 
experimental mean squared estimation errors in T40, T50, Fn 
and SmLPC are reduced 76%, 82%, 80% and 63%, respec- 
tively. It can also be observed that the SVD tuner selection 
approach, which is designed to reduce the estimation error 
55 bias, does in fact reduce the sum of squared biases relative to 
the subset of health parameters approach. However, the SVD 
approach is also found to increase the estimation variance, 
which contributes to its overall mean squared estimation 
error. 

60 A visual illustration of the effect that tuner selection has on 
Kalman filter estimation accuracy can be seen in FIGS. 6-9, 
which show actual and estimated results for the auxiliary 
parameters T40, T50, Fn and SmLPC respectively. Each plot 
shows a 300 s segment of the evaluated test cases. The step 
65 changes that can be observed in each plot every 30 s corre- 
spond to a transition to a different health parameter vector. 
True model auxiliary parameter outputs are shown in black, 
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and Kalman filter estimates are shown in red. In each figure 
the information is arranged top to bottom according to tuner 
selection based upon: a) a subset of health parameters; b) 
singular value decomposition; and c) the new systematic 
selection strategy. The information shown in these figures 
corroborates the information in FIG. 5; namely all three tuner 
selection approaches produce unbiased estimates of T40 
(FIG. 6), while the systematic tuner selection strategy yields 
a noticeable reduction in the total squared estimation error 
(squared bias plus variance) of all four auxiliary parameters. 

The presented systematic tuner selection strategy mini- 
mizes the mean squared error of the on-line estimator at 
steady-state operating conditions, taking advantage of prior 
knowledge of engine health parameter distributions. As such 
it is somewhat analogous to the maximum a posteriori (MAP) 
estimation method which is commonly applied for ground- 
based aircraft gas turbine engine gas path analysis. This leads 
to the question, how does the on-line Kalman filter estimation 
accuracy compare to MAP estimation accuracy? Prior to 
making this comparison the mathematical formulation of the 
MAP estimator is briefly introduced. Here a steady-state 
model of the Measurement process in the following form is 
applied 

y k =Hh k +v k ( 39 ) 

where the matrix H relates the effects of the health parameter 
vector, h, to the sensed measurements, y. From Eq. (19), it can 
be seen that H is equivalent to C(I-A) -1 L+M. The maximum 
a posteriori (MAP) estimator follows the closed Form expres- 
sion 

k k =(P h - 1 +H T R- l H)- 1 H r R~ l y k ( 40 ) 

The MAP estimator is capable of estimating more 
unknowns than available measurements due to the inclusion 
of a priori knowledge of the estimated parameter covariance, 
P /r However, the MAP estimator, unlike a Kalman filter, is not 
a recursive estimator and does not take advantage of past 
measurements to enhance its estimate at the current time step. 
Furthermore, the MAP estimator only considers a static rela- 
tionship between system state variables and measured out- 
puts — it does not consider system dynamics. Because of these 
differences a Kalman estimator with optimally selected tun- 
ing parameters should outperform the MAP estimator How- 
ever, under steady- state conditions, with minimal sensor 
noise the two estimation approaches should produce similar 
results. To test this theory, a MAP estimator was designed and 
its estimation accuracy was compared to a Kalman filter with 
tuning parameters optimally selected to minimize the estima- 
tion errors in the health parameter vector h. First, the two 
estimators were designed and evaluated using the original 
sensor noise levels shown in FIG. 3. Next, the sensor noise 
levels were set to V^oth of their original levels, the estimators 
were re-designed, and the comparison was repeated. Monte 
Carlo simulation evaluations as previously described were 
applied (i.e., 375 random health parameter vectors, 30 s in 
duration, with estimation accuracy calculations based upon 
the last 10 s of each 30 s test case). Theoretical and experi- 
mental estimation errors are shown in FIGS. 10 and 11 for the 
original noise and reduced noise levels, respectively. At origi- 
nal noise levels the Kalman estimator is able to produce 
smaller estimation errors. However, at the reduced noise level 
the two estimation approaches are found to be nearly identi- 
cal. This comparison validates that the Kalman estimation 
approach is indeed producing a minimum mean squared esti- 
mation error as intended, while providing the capability to 
support real-time on-line estimation under dynamic operating 
scenarios. 
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While the systematic tuner selection approach presented 
here appears promising for on-line Kalman filter based 
parameter estimation applications, there are several practical 
considerations which need to be assessed when applying such 
5 a technique. The optimization routine attempts to minimize 
the overall squared estimation error — both bias and vari- 
ance — under steady-state operating conditions. The minimi- 
zation of the estimation variance in particular can come at the 
expense of dynamic responsiveness of the Kalman filter. To 
to illustrate this consider the time history plots of actual versus 
estimated T40 shown in FIG. 12. The top plot shows Kalman 
filter estimation results using a tuning parameter vector sys- 
tematically selected to minimize the error in four auxiliary 
parameters (T40, T50, Fn and SmLPC) as presented in the 
15 previous section. The bottom plot shows Kalman filter esti- 
mation results using a tuning parameter vector systematically 
selected to minimize the estimation error in T40 only. At time 
100 s a step change in the health parameter input vector is 
introduced into the engine model; this allows the dynamic 
20 response of the two estimators to be compared. It can be 
observed that T40 estimation variance in the bottom plot is 
reduced, as is the mean steady- state estimation error (>300 s). 
This is not surprising since one would generally expect 
improved results when optimizing to minimize the error in a 
25 single parameter, as opposed to multiple parameters. How- 
ever, the estimator shown in the bottom plot does require a 
significantly longer time to reach steady-state convergence. 
Conversely, the estimator designed to minimize the steady- 
state error in four auxiliary parameters (top plot) is unable to 
30 place as much emphasis on T40 estimation variance reduc- 
tion, but it is able to track dynamic changes in T40 more 
rapidly. This example illustrates the inter-dependence 
between estimation variance and responsiveness. Therefore, 
it is prudent for a designer to evaluate the Kalman filter to 
35 ensure that it tracks engine dynamics acceptably. If the 
dynamic response is unacceptable, the optimization routine 
can be re-run placing more weight on estimation error bias 
reduction, and less weight on variance reduction. 

In an embodiment, the present approach produces an opti- 
40 mal set of tuning parameters, not just at a single operating 
point but rather a globally optimal tuning parameter vector 
universally applicable over the range of operating conditions 
that an engine is expected to experience. A potential approach 
to selecting a single “globally optimal” tuning parameter 
45 vector is to modify the optimization routine to minimize the 
combined estimation error over multiple engine operating 
points such as takeoff, climb and cruise. This would be a 
straightforward modification to the Matlab® optimization 
routine, but it would increase the computational time required 
50 to calculate the result. Since the systematic tuner selection 
process is only envisioned to be done once during the system 
design process, this will not impact the on-line execution 
speed of the Kalman filter. It is anticipated that the application 
of globally optimal tuners will result in some estimation 
55 accuracy degradation relative to tuners optimized for indi- 
vidual operating points, although this has not yet been veri- 
fied or quantified. 

A systematic approach to tuning parameter selection for 
on-line Kalman filter based parameter estimation has been 
60 presented. This technique is specifically applicable for the 
underdetermined aircraft engine parameter estimation case 
where there are fewer sensor measurements than unknown 
health parameters which will impact engine outputs. It creates 
and applies a linear transformation matrix, V*, to select a 
65 vector of tuning parameters which are a linear combination of 
all health parameters . The tuning parameter vector is selected 
to be of low-enough dimension to be estimated, while mini- 
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mizing the mean-squared error of Kalman filter estimates. 
The multiparameter iterative search routine applied to opti- 
mally select V* was presented. Results have shown that while 
the transformation matrix returned by the optimization rou- 
tine is not unique (different matrices can be found which 
produce a global minimum of the objective function), the 
routine is effective in returning a transformation matrix which 
is optimal, or near optimal, regardless of its initial starting 
guess of the matrix. The efficacy of the systematic approach 
to tuning parameter selection was demonstrated by applying 
it to parameter estimation in an aircraft turbofan engine linear 
point model. It was found to significantly reduce mean 
squared estimation errors compared to the conventional 
approach of selecting a subset of health parameters to serve as 
tuners. In some parameters the mean squared estimation error 
reduction was found to be over 80%. These estimation 
improvements were theoretically predicted and experimen- 
tally validated through Monte Carlo simulation studies. 

The systematic approach to Kalman filter design may be 
applicable for a broad range of on-board aircraft engine 
model-based applications which produce estimates of 
unmeasured parameters. This includes model based controls, 
model -based diagnostics, and on-board life usage algorithms. 
It also may have benefits for sensor selection during the 
engine design process, specifically for assessing the perfor- 
mance estimation accuracy benefits of different candidate 
sensor suites. 

The methods and processes described herein may be per- 
formed by a hardware or software system or combination 
thereof. For example, a computing system may be configured 
to receive input signals from an engine by way of a plurality 
of sensors, and model the health parameters using the Kalman 
filter approach described herein. The computing system may 
comprise a processor or CPU, a memory, and a storage 
device, as is known in the art. The computing system may 
further include a logic configured to receive the input signals 
from the engine and determine estimates for unknown param- 
eters. The computing system may modify the engine control 
based on an estimated unknown parameter or parameters 
determined by the logic using the methods and processes 
described herein. 

As used herein, the term “logic” includes but is not limited 
to a software, a firmware, an executable program, a hardware 
or hard-wired circuit, or combinations thereof. For example, 
based on a desired application or needs, logic may include a 
software controlled microprocessor, discrete logic like an 
application specific integrated circuit (ASIC), an analog cir- 
cuit, a digital circuit, a programmed logic device, a memory 
device containing instructions, or the like. Logic may include 
one or more gates, combinations of gates, or other circuit 
components. Logic may also be fully embodied as software. 
Where multiple logical logics are described, it may be pos- 
sible to incorporate the multiple logical logics into one physi- 
cal logic. Similarly, where a single logical logic is described, 
it may be possible to distribute that single logical logic 
between multiple physical logics. 

Although the preferred embodiments of the present inven- 
tion have been illustrated in the accompanying drawings and 
described in the foregoing detailed description, it is to be 
understood that the present invention is not to be limited to 
just the preferred embodiment disclosed, but that the inven- 
tion described herein is capable of numerous rearrangements, 
modifications and substitutions without departing from the 
scope of the claims hereafter. 

We claim: 

1. A method of optimizing the estimated performance and 
estimated condition of an engine comprising: 
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estimating a plurality of engine parameters based on a 
plurality of sensor inputs, wherein said engine param- 
eters include a set of health parameters and a set of 
performance parameters, and wherein the number of 

5 health parameters in said set is greater than the number 
of sensors; 

determining a relationship between said set of health 
parameters and a tuning vector, wherein said tuning 
vector is a linear combination of said entire set of health 
parameters; 

determining at least one state equation in terms of said 
tuning vector; 

solving said at least one state equation to determine an 

1 5 estimate for said set of health parameters ; 

estimating said performance parameters based on said 
health parameter estimates; 

calculating the error in said estimation of said health 
parameters and said performance parameters; 

20 optimizing said relationship between said tuning vector 
and said set of health parameters based on said error 
calculation; 

determining an optimized tuning vector based on said opti- 
mized relationship; 

25 determining at least one optimized state equation in terms 
of said optimized tuning vector; 

solving said at least one optimized state equation to deter- 
mine an optimized estimate of said engine parameters. 

2. The method of claim 1, wherein calculating the error 

30 includes calculating the mean estimation error. 

3. The method of claim 1, wherein said performance 
parameters are a function of said health parameters. 

4. The method of claim 1, wherein said relationship 
between said tuning vector and said set of health parameters 

35 is optimized to minimize the error in estimating a given sub set 
of said health parameters. 

5. The method of claim 1, wherein said relationship 
between said tuning vector and said set of health parameters 
is optimized to minimize the error in estimating a given subset 

40 of said performance parameters. 

6. The method of claim 1, wherein said relationship 
between said tuning vector and said set of health parameters 
is optimized to minimize the error in estimating a combina- 
tion of a given subset of said health parameters and a given 

45 subset of said performance parameters. 

7. The method of claim 1, wherein the dimension of said 
tuning vector is less than or equal to the number of said 
sensors. 

8. The method of claim 1, wherein said relationship 

50 between said health parameters and said tuning vector is 

defined by a matrix. 

9 . The method of claim 8, wherein a first dimensions of said 
matrix is equal to the dimension of said tuning vector. 

10. The method of claim 9, wherein a second dimension of 

55 said matrix is less than or equal to the number of said health 

parameters. 

11. The method of claim 1 further comprising the step of 
modifying control of said engine based on said optimized 
estimate of said engine parameters. 

60 12 . The method of claim 1 wherein solving said state equa- 

tion includes constructing a Kalman filter. 

13. An optimal tuning system comprising: 

a plurality of sensors capable of sensing engine param- 
eters; 

65 a logic in operable communication with said plurality of 
sensors and configured to receive input from said plu- 
rality of sensors, said logic configured to determine opti- 
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mized parameters related to the estimated performance 
and condition of said engine; 
wherein said logic estimates a plurality of engine param- 
eters based on said sensors, said engine parameters 
including a set of health parameters and a set of perfor- 
mance parameters, and wherein the number of health 
parameters in said set is greater than the number of 
sensors; 

wherein said logic determines a relationship between said 
set of health parameters and a tuning vector that is a 
linear combination of said entire set of health param- 
eters; 

wherein said logic determines at least one state equation in 
terms of said tuning vector and solves said at least one 
state equation to determine an estimate for said set of 
health parameters; 

wherein said logic estimates said performance parameters 
based on said health parameter estimates and calculates 
the error in said estimation of said health parameters and 
said performance parameters; 
wherein said logic optimizes said relationship between 
said tuning vector and said set of health parameters 
based on said error calculation; 
wherein said logic determines an optimized tuning vector 
based on said optimized relationship; 
wherein said logic determines at least one optimized state 
equation in terms of said optimized tuning vector 
wherein said logic solves said at least one optimized state 
equation to determine an optimized estimate of said 
engine parameters. 
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14. The optimal tuning system of claim 13, wherein said 
logic comprises multiple logics. 

15. The optimal tuning system of claim 14, wherein said 
optimized tuning vector is defined within one of said multiple 

5 logics. 

16. The optimal tuning system of claim 13, wherein said 
logic calculates the mean estimation error in said estimation 
of said health parameters and said performance parameters. 

17. The optimal tuning system of claim 13, wherein said 

10 performance parameters are a function of said health param- 
eters. 

18. The optimal tuning system of claim 13, wherein said 
logic optimizes said relationship between said tuning vector 
and said set of health parameters to minimize the error in 

15 estimating a given subset of said health parameters. 

19. The optimal tuning system of claim 13, wherein said 
logic optimizes said relationship between said tuning vector 
and said set of health parameters to minimize the error in 
estimating a combination of a given subset of said health 

20 parameters and a given subset of said performance param- 
eters. 

20. The optimal tuning system of claim 13, wherein said 
logic solves said state equation by constructing a Kalman 
filter. 

25 21. The optimal tuning system of claim 13, wherein said 

logic is capable of modifying control of said engine based on 
said optimized estimate of said engine parameters. 

22. The optimal tuning system of claim 13, wherein the 
dimension of said tuning vector is less than or equal to the 

30 number of said sensors. 



